The Holographic Superconductor Vortex 
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A gravity dual of a superconductor at finite temperature has been recently proposed. We present 
the vortex configuration of this model and study its properties. In particular, we calculate the 
free energy as a function of an external magnetic field, the magnetization and the superconducting 
density. We also find the two critical magnetic fields that define the region in which the vortex 
configurations are energetically favorable. 
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INTRODUCTION 

The Gauge/Gravity duality, that relates strongly in- 
teracting gauge theories to theories of gravity in higher 
dimensions, has opened a new window to study many dif- 
ferent strongly interacting systems. The applicability of 
this approach is very vast ranging from particle physics 
to plasma and nuclear physics. In Ref. [1| a model for a 
dual description of a superconductor was proposed. The 
model showed to have a critical temperature under 
which the system goes into a superconducting phase. The 
properties of this phase have been thoroughly studied , 
showing a resemblance with those of a Type II supercon- 
ductor. In spite of this, Abrikosov vortices, known to 
happen in Type II superconductors, have not yet been 
obtained. The purpose of this letter is to show that in 
this type of gravity duals vortex solutions indeed exist 
and can be energetically favorable in the presence of ex- 
ternal magnetic fields. Due to the nonlinear nature of 
these configurations, we will have to rely on numerical 
methods. Among other physical properties, we will cal- 
culate the free energy and the range of the magnetic field 
Sci < S < i?c 2 at which the superconductor is at the 
intermediate phase (Shubnikov phase) characterized by 
vortex configurations. Further aspects of these solutions 
will be presented elsewhere. 



THE MODEL 

The physical system to study is a conformal strongly 
coupled superconductor in 3D at finite temperature and 
charge density. Its gravitational dual theory [ij is an 
asymptotically AdS-Schwarzschild space-time in 4D. The 
gravitational degrees of freedom are coupled to an U(l) 
gauge field and a complex scalar ^ . The action that 
summarizes the above model is given by 




5* 



with C^-F'^ 
4 



_J_(i? + A)-l£ 



Gat is the 4D gravitational Newton constant, the cosmo- 
logical constant A defines the asymptotic AdS radius L 
via the relation A = — 3/L^ and — — iA^. We 
use the convention where the metric G has signature 
(— , +, +, +), with coordinates (t, z, r, 0) where t is time, z 
is the holographic direction such that the AdS-boundary 
occurs at ^ = 0, and (r, (j)) are polar coordinates parame- 
terizing the remaining 2D plane. For the scalar mass 
we will focus on two possible values: = — 2, 0. Other 
values are expected to give similar behaviors [3|. 

We will work in the so-called probe approximation, 
where the gravity sector is effectively decoupled from the 
matter sector and therefore, there is no back-reaction on 
the background metric due to C. This regime is achieved 
in the limit of large when compared to the gravitational 
strength. In this limit we can, without loss of generality, 
fix g = 1. In our conventions, the background AdS- 
Schwarzschild Black hole (BH) metric is given by 
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— {-f{z)dt' 



dr^ 



z^fiz) 



dz^ 
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where J{z) = 1 - {z/zhY- 

As we are considering the theory at finite temperature, 
we have to take the Euclidian regime with compact time 
it e [0, 1/T] where T = S/iAnzh). Therefore, the holo- 
graphic coordinate runs from the AdS-boundary at ^ = 
to the BH horizon at z = . Notice that we work with a 
planar BH with energy per unit area e = / {%ttG m zf^) ■ 
Then, the AdS/CFT duality tells us that the above are 
precisely the temperature and energy density of the dual 
superconductor. 

The gauge field has the usual AdS-boundary behavior 
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where a,y = (/i, a^) corresponds to the potentials on the 
dual CFT, while J^, = (— p, Jj) plays the role of the con- 
jugated currents. We will consider the case in which the 
charge density p is fixed constant. The other potentials 
ai are related to turning on either electromagnetic fields 
or sample velocities in the dual CFT, depending on the 
interpretation we give to the AdS / CFT duality. The first 
interpretation is what we will use in this article, while the 
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second one is relevant for superfluids ^ Similarly, the 
scalar field has the following AdS-boundary behavior 
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where A = 2, 3 (for — —2, 0) corresponds to the di- 
mension of the dual operator Oa responsible for the U(l) 
breaking, and b determines the vacuum expectation value 
of this operator. The value of a corresponds to an explicit 
breaking of the U(l) symmetry and will then be turned 
to zero ^ . Having fixed , the only parameters of the 
model are the scales T and y/p. 

It has been reported in Ref. [J, that for p 7^ the 
system undergoes a phase transition at 



Tc ~ 0.12^ for TO^ = -2, 
Tc ~ 0.09VP for TO^ = , 



(5) 



where the two phases are related to a charged BH and a 
charged BH with a non-trivial scalar hair. At T < Tc, the 
system is at the hairy phase corresponding to a supercon- 
ducting phase. In Refs. the model was also studied 
in the presence of an external magnetic field B using a 
dyonic BH with a probe scalar field. The result was a 
bounded superconducting region or drop, that squeezes 
to zero size as we increase B. The above suggested that 
we are dealing with a Type II superconductor. If this 
is the case, Abrikosov vortex configurations should be 
present in this model. 

We stress that, as is usual in this approach, we are 
treating the electromagnetic field of the 3D dual theory 
as a nondynamical background. This corresponds to take 
the 3D electric charge e — > 0, while keeping constant B 
and p. 



THE VORTEX SOLUTION 



We use the Ansatz given by 

M' = V(r,^)e»^ Ao^Ao{r,z), ^ A^{r, z) , (6) 

with all other fields set to zero. This Ansatz preserves 
global U(l) transformations when combined with a ro- 
tation in the 2D plane. The fields Ar,Az can be con- 
sistently set to zero since our Ansatz fulfills drArg['^] = 
dzArg[^] = 0. The winding number n G Z determines 
different topological solutions. With the above Ansatz we 



^ In fact, the vortex solution wc present in this article can be 

identified with vortex configurations in a superfluid, once the 

appropriated reinterpretations are made. 
^ For the case = —2 there is the possibility to have 6 = and 

a ^ corresponding to have a dual CFT operator of dimension 
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obtain from Eq. ((T|) the following equations of motion: 
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d. ifdzA^) +rdr[ -drA^) - ^(A^ - n) = , 
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2V'' 



fdiAo + -Or (rdrAo) - -V^o = 



(7) 



In order to describe a dual superconductor at fixed p 
in the presence of an external magnetic field B, the 
AdS/CFT correspondence tells us that we must impose 
the AdS-boundary conditions 

V'U=o = , dzAoU=o = -p , A^U=o = \r^B , (8) 

for the case m? = 0, while for m? — —2 the first condition 
must be 9zV'|z=o = (this is equivalent to set a = in 
Eq. (HI)). At the horizon z = we require the field con- 
figurations to be regular; in particular we set Af)\z=Zh = 
as usual, to have a well-defined Euclidean continuation. 
Similar reasoning at r = implies that for rt ^ 



0, 



(9) 



while for n = 0, 9r''0|r=o = 0. We will be considering a 
3D superconductor of radius R that we will take to be 
much bigger than the vortex radius. This is implemented 
by setting a nonzero p extending from r = to r = i?. 

The 2D system of the three partial differential equa- 
tions of Eq. ^ is nonlinear, and therefore requires to be 
solved numerically. For this purpose we have used the 
COMSOL 3.4 package 0|- In our numerical studies we 
have chosen 
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This corresponds to 



T 



T = 0.065 Vp ■ 



0.74 (0.54), 



(10) 



(11) 



for the case of = (—2). 

In Fig. 1 we show the order parameter (Oa) — 
\z^~^dz^\z=o of the dual superconductor. We can see 
that this goes to zero at the origin where the vortex is 
placed. For the value of the magnetic field, we have cho- 
sen 

2n 

i?2 ' 



Bn 
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corresponding to the value at which the magnetic flux 
crossing a surface of constant z, $ = /d0 J^^rdrB, 
equals 2nn. This is the quantized flux going through 
the n-vortex of the dual superconductor. 
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FIG. 1: Order parameter (Oa) for the n = 1 (solid) and 
n — 2 (dashed) vortex configuration. The lower (upper) curves 
correspond to the case — (—2). Presented in units of 



FREE ENERGY, MAGNETIZATION AND 
CRITICAL MAGNETIC FIELDS 



We are interested to deterniine the free energy of the 
superconductor configurations with n = 0,1,2 to know 
which one is energetically favorable as we vary B. By 
the AdS/CFT, the free energy F of the superconductor 
is given by 



T 



= Sp 



TT 

T 



drrAodzAo 



(13) 



z=0 



where the right-hand side is evaluated on-shell in the 4D 
theory with the boundary conditions given in Eq. ([8]). 
The second term of Eq. (fT3|) has been added to guaran- 
tee the variational principle when working at fixed dzAo 
on the AdS-boundary. Since, as we will see, the phase 
transition to vortex configurations occurs at small values 
of -B, we can treat the magnetic field as a small pertur- 
bation and separate the solution as 

TA^V^ + fJVj, Ao^Ao + 6Ao, A^~^ A^ + SA^, (U) 

where the unperturbed solution (jp, Aq, A^) corresponds 
to that at zero external magnetic field, i.e., A^j^^o = 0, 
while the perturbation {Sip^ SAq, SAfj,) must fulfill 

6A^U=o = ^r^B , 9.MoU=o = , 5VU=o = , (15) 

for = and dz5il^\z=o = for = —2. By integrat- 
ing by parts the free energy of the n- vortex configuration 
can be written, up to terms, as 



F„(B) ~ F„(0) - a„B + -/3„B2 ^ 



(16) 



where we have defined 

F„(0) = 27r / dr 

Jo Jo 

— TT drrAo OzAq 
Jo 

2n f^dr^^ „ ^ 
ck„ = — / — dAA,o,A 



z=0 



B 



2tt dr 
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SAsdJA 



(17) 



Notice that the positive-defined quantities a„ and /?„ do 
not depend on B, since SA^ cx SA^\z=o oc B. Eq. (flB)) 
has a simple interpretation in terms of the magnetization 
M of the superconductor. Using M = —dF/dB, we can 
write 



Fn{B)=F„{0)- [ M„dB, 
Jo 



(18) 



where the magnetization of the n-vortex configuration 
Mn in the z-component is given by 



idrr{rx J)z — drrj, 



From the AdS/CFT dictionary, we have that 
SF 



SA'P\z^o 



= dzA^ + dJA4 



(19) 



(20) 



z=0 



that together with Eq. leads to our final expression 
for the magnetization 



(21) 



Using this expression into Eq. p^ . we recover the free 
energy of Eq. pB]) . 

For the free energy at i? = we obtain 

FniO)^Fo{0) + 0.9il.5)n'ln[Rp'^'']^ + c„, (22) 

where cq — 0, ci ~ 1.2(3.7)y^, C2 — 0.3(4)y/p and 

Fo(0)~5(4)i?V7^, (23) 



for the case 



0(— 2). This shows that, as expected, 



the vortex configurations have for B = a larger energy 
than the n — solution. Note that Fo{0) grows with the 
volume of the superconductor (oc R^), although not the 
difference i^i,2(0) — Fo(0) that is only logarithmically sen- 
sitive to Rior R ^ oo, as expected for 3D vortices in the 
absence of electromagnetic fields. For the magnetization 
we find 

a„ ~ 0.4(0.7) , /3„ ~ 0.05(0.09)i?'*VP • (24) 
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From Eq. ((16)) it is clear that there is a critical value 
for B at which the difference between the free energies 
Fi{B) — Fq{B) is zero. This value is usually referred as 
Be I and marks the beginning of the mixed phase where 
the magnetic field starts to penetrate the superconductor. 
For the case of = we have 



fi(g)-fo(g) 

VP 



0.91n[i?p 



1/21 



1.2-0.8^, (25) 



that for R — 5Q/yfp equals to zero at 



B. 



c 1 



6Bi 



(26) 



where Bi is defined in Eq. p^ . For = — 2 we get 
similar values, Bd — 7Bi. At higher magnetic field val- 
ues than Bel the vortex configuration is preferred. No- 
tice that for i? ^ oo, we have Bi and therefore 
Bel — > 0, indicating that the non- vortex solution is never 
favorable at any B ^ 0. 

For the configuration with n = 2, we find that its free 
energy is less than that for n = 0, 1 if B > 10(14)_Bi 
for — 0(— 2). At this high magnetic field, however, 
we expect that the free energy of a solution with two 
n — I vortices will be energetically more favorable, as 
it happens in Type II superconductors. Indeed, for two 
vortices sufficiently separated we expect 



F{B) ~ i;^o(0) + 2[(^^i(0)-^^o(0))-aiB] 
+ Eint + lpiB\ 



(27) 



where -Eint is the interaction energy between the two vor- 
tices. Therefore the difference between the free energy 
of two n = I vortices and one n — 2 vortex goes as 
AF ~ Eint - 1.8(3) ln[i?pi/2]^ for = 0(-2). As 
a consequence a configuration with two n = I vortices 
will be preferred for Ei^t < 1.8(3) ln[i?pi/2]^ ^j^a^ 
expected for a large superconductor. 

On the other hand, as B increases from Bci, a con- 
figuration with more and more vortices is expected to 
be favorable, until we reach a certain critical value i?c2 
at which there is another phase transition; for B > Bc2 
the normal phase is preferred. We estimate this value by 
the magnetic field at which the superconducting region 
of the n = 0,1 configurations shrink to zero size. We find 
Be2 ^ 3(5)p for m2 = 0(-2). 

In Fig. 2 we plot the values of the free energy as a 
function of B for the configurations n = 0, 1, 2 from the 
exact numerical solutions. We can see that the critical 
magnetic values at which the lines cross are similar to 
the approximate ones given above. 

Finally, we calculate the "superconducting density" 
ns{r) defined as 



SF 



(28) 
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FIG. 2: Free energy for the m? = case as a function of the 
external magnetic field for the n = (solid), n = 1 (dashed) 
and n = 2 (dotted) vortex configuration. Presented in units 
of^=l- 



where in the last equality we have used Eq. In 
Fig. 3 we show risir) for the different configurations. 
We notice that the vortex configuration fuUfiUs (J^) — 
—ns{r){5A^\z=o — n), as expected from a spontaneously 
broken U(l) symmetry. For a non- vortex configura- 
tion the superconducting density is constant ns(r) ~ 
0.28(0.48)^ for = 0(-2). This determines the pene- 
tration length X = 1/ {e^/ni) where e is the electric charge 
of the dual superconductor. 
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FIG. 3: Superconducting density na{r) for the n = 1 (solid) 
and n = 2 (dashed) vortex configuration. The lower (upper) 
curves correspond to the case = (—2). Presented m 
units of = 1 . 

Note Added: While finishing this paper, we learned 
of Ref. which has also studied the vortex solution in 
holographic superconductors. 
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